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In this paper we present an analysis to determine the existence of singularities in spatially ho- 
mogeneous anisotropic universes filled with nonlinear electromagnetic radiation. These spaces are 
conformal to Bianchi spaces admitting a three parameter group of motions G3. For these models 
we study geodesic completeness. It is shown that with nonlinear electromagnetic field some of the 
Bianchi spaces are geodesically complete, like G3IX and G3VIII; however, completeness depends on 
the curvature of the space. When certain topology is assumed, Bianchi G3IX presents geodesies 
that are imprisoned. It is surprising that in the linear limit (Maxwell field) the spacetimes are 
singularity-free even if the curvature parameter is zero. 
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I. INTRODUCTION 

General Relativity breaks down at certains points 
called singularities. The singularity theorems, that were 
shown for the first time 30 years ago by Stephen Hawk- 
ing, Robert Geroch and Roger Penrose Q , gave the bad 
news that under certain conditions the singularities of 
the spacetime would appear in the future or in the past. 
An important subject that is not yet determined is the 
nature of the singularities of the spacetime in general 
relativity. 

There are several criteria that tell us if a spacetime is 
singular, such as energy conditions or the divergence of 
the curvature invariants. However, these conditions are 
only necessary but not sufficient. If one of them fails to 
hold, one is ignorant if the spacetime possesses a singular- 
ity or not. The common characteristic of singular space- 
times is the existence of incomplete causal geodesies. We 
shall adopt as definition of the absence of singularities 
the completeness of causal geodesies (g-completeness) , i. 
c. no observer in free fall leaves the spacetime in a fi- 
nite proper time. We note, however, that not even the 
geodesic completeness guarantees the absence of singular- 
ities (see a counterexample by Geroch 0, of a geodesi- 
cally complete spacetime with incomplete non-geodesic 
curves) . 

Mars and Senovilla presented a Bianchi IX model with 
an energy- momentum tensor of perfect fluid type 3] that 
is completely regular everywhere in the energy density 
and presure and is geodesically complete. This model 
provides an example of how the presence of matter can 
help to avoid the initial singularity while keeping physical 
compelling dominant energy conditions. 

Recently, more studies about singularity-free Bianchi 
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type IX models have arisen; for example, Berger used the 
method of consistent potentials Q to explain how min- 
imally coupled (classical) scalar field can suppress mix- 
master oscilations [j| when approaching the singularity of 
Bianchi IX cosmology. Toporensky and Ustiansky stud- 
ied the possibility of non-singular transition from con- 
traction to expansion in the dynamics of the Bianchi IX 
cosmological model with minimally coupled massive real 
scalar field Q. 

In this paper we study singularities in spatially ho- 
mogeneous but anisotropic spaces where the energy- 
momentum tensor corresponds to (source-free) nonlinear 
electromagnetic field, in particular, Born-Infeld field. For 
these spaces the energy momentum tensor satisfies the 
dominant energy condition and the matter does not di- 
verge; three cases possess regular Weyl invariants. These 
conditions point to the absence of singularities provided 
the spaces do pass the test of geodesic completeness. The 
asumed spacetime includes particular cases of Bianchi 
spaces. The analysis of the geodesies lead us to con- 
clude that some of the studied Bianchi spaces avoid the 
singularity, like Bianchi G3IX and G3VIII, while Bianchi 
G3II might present singularity. Spaces like Bianchi G3III 
and G3I as well as a particular case of Kantowski-Sachs 
present singularity in the curvature, i.e. the invariants 
diverge after a finite time, however, in two cases the 
geodesies are complete. 

The solutions investigated in this papers have already 
been discussed. Approximate solution was presented in 
7] and the posibility of (NLED) as source of inflation 
was addressed in |8j| . The singularity structure is studied 
in the present paper. 

Before getting into the main subject, we give a sum- 
mary of non-linear electrodynamics (NLED) in Sec. II. 
In Sec. Ill the studied solutions of the coupled Einstein- 
Born-Infeld (EBI) equations are presented along with 
comments on their geometric structure, energy condi- 
tions and scalar invariants. In Sec. IV we analyze the 
geodesies and the nature of the singularities that appear 
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in the models. In Sec. V we address the limit of lin- 
ear electrodynamics (Einstein-Maxwell) and it turns out 
that for this field all solutions are geodesically complete. 
In the final section we draw some conclusions. 



II. BORN-INFELD NONLINEAR 
ELECTRODYNAMICS 

One of the motivations to explore spacetimes with non- 
linear electromagnetic fields is that in early epochs of 
the universe the magnetic fields exceeded some 10 15 G. 
At such values of the electromagnetic field, the interac- 
tion of photons with themselves becomes important and 
classical electrodynamics is not valid anymore [jj. To 
overcome the use of quantum electrodynamics (QED) we 
shall employ the Born-Infeld theory. 

The nonlinear electrodynamic theory presented by 
Born and Infeld (1934) possesses a Lagrangian which 
takes into account quantum electrodynamical correc- 
tions, this fact was noted by Euler and Heisenberg soon 
after the proposal by Born and Infeld and later on it was 
remarked by Schwinger. In the papers by Heisenberg 
and Euler [l(| and by Drummond and Hathrell it 
was demonstrated that the Born-Infeld lagrangian has 
the same dependence in the electromagnetic invariants 
than QED for the one-loop approach. On the other side, 
based on very general considerations, Schwinger |l2j | ob- 
tained an expression for the lagrangian that coincides, 
within numerical factors of order unity, with the BI la- 
grangian in the same approximation. Therefore, EBI the- 
ory should give us, at least qualitatively, characteristics 
of early universes that could not be obtained within clas- 
sical electrodynamics. 

Moreover, lately the BI type action has arisen in 
string theory: it happens that open strings with Dirichlct 
boundary conditions (D-branes) are described by an ef- 
fective action of Dirac-Born-Infeld, that is related to the 
effective action of BI type in the theory of open strings. 
The action to low energies corresponding to these objects, 
called branes, can be written like a term of BI plus a term 
of the Weiss-Zumino or Chern-Simmons type jl3]. The 
BI strings are solitons that represent terminal strings on 
branes. 

On the other side, in relation with cosmological models 
with electromagnetic fields, Brill studied Einstein- 
Maxwell solutions in a homogeneous and nonisotropic 
universe. His solution is a generalization of Taub's vac- 
uum universe and represents a closed universe of topology 
Rx S 3 filled with gravitational and electromagnetic radia- 
tion. This universe expands in two spatial directions and 
contracts along the third one; moreover its invariants are 
finite. Also Hughston and Jacobs 0] investigated the 
behavior of source-free, homogeneous, electromagnetic 
fields in Bianchi type cosmologies. They showed that 
there can be no pure electric (magnetic) fields in types 
VIII, IX, V, VI (h ^ 0) and VI (h ^ -1). Types III, VI 
(h = — 1, h = Q) allow one independent component of the 



magnetic field. Type II allows two independent magnetic 
field components while type I puts no constraints on the 
magnetic field components. Note that results in were 
obtained for Maxwell fields in the Bianchi geometry, non 
coupled to Einstein equations. In a previous paper, Ha- 
cyan 0] investigated the effect of BI electrodynamics in 
the initial singularity of a Bianchi I space. It turned out 
in that case that the presence of a BI field smoothed the 
oscillations of the metric functions when approaching the 
singularity. In the cases addressed in the present paper, 
Bianchi types G3I, G3III as well as Kantowski-Sachs ad- 
mit an electric or magnetic non linear field (not both) and 
G3II, G3VIII and G3IX admit both components, more- 
over, the presence of the electromagnetic field prevents 
the occurrence of singularities in some Bianchi spaces. 

A. The Born-Infeld field parameterization 

We describe in what follows the main characteristics 
of the BI formalism used to determine the EBI solutions. 

The Born-Infeld (BI) nonlinear electromagnetic theory 
is self-consistent and satisfies all natural requirements. 
Its lagrangian depends on the BI constant b, the max- 
imum field strenght; it also depends in nonlinear form 
on the electromagnetic invariants, P and Q, of the the 
antisymmetric tensor P^ v , the generalization of the elec- 
tromagnetic tensor F^ v ; the BI lagrangian is given by 

L BI = -±P'"'F flv + H(P,Q) (1) 

where Ti.(P, Q) is the so called structural function. The 
structural function is constrained to satisfy some physi- 
cal requirementes: (i) the correspondence with the linear 
theory of Maxwell in the limit when b goes to infinity 
(H(P,Q) = P + 0(P 2 ,Q 2 )); (ii) the parity conserva- 
tion {H{P,Q) = H(P,-Q)); (hi) the positive definite- 
ness of the energy density (Tt,p > 0) and, (iv) the re- 
quirement of the timclikc nature of the energy flux vec- 
tor {PH, p +QH,q-H < 0). Conditions (hi) and (iv) 
amount to the fulfilment of the dominant energy condi- 
tion (DEC) Q: T^yU^vy > 0, for every time direction 
it M , (u^u^ < 0) and T^u^ is a nonspacelike vector. The 
energy-momentum tensor in terms of the structural func- 
tion is given by: 

= -U,p P^ a P a v + g MV (2PH lP +QH, Q -H) . (2) 

The structural function for the BI field is given by 

H(P, Q) = b 2 - Vb 4 -2b 2 P + Q 2 , (3) 

In the linear limit, which is obtained by taking b — * 00, 
then Ji = P and P^ v — F^. P^ u and F^ are related 
through the material or constitutive equations: 

F^=H,pP^+H, Q P^, (4) 
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where P^ — —je.^af}P a ^ is the dual of P M „. An exten- 
sive treatment on nonlinear electrodynamics in curved 
spacetimes can be consulted in Il7l . 

In the null tetrad formalism |l8|. and dealing with a 
Petrov type-D metric, one can always align the direc- 
tions of the real null vectors e 3 and e 4 along the double 
Debever-Penrose (DP) vectors. We shall assume that the 
eigenvectors of F a b (and consequently the ones of P a b) are 
also aligned in the direction of the DP vectors. Hence the 
novanishing components of F a f, are F\2 and -F34 (P12 and 
^34) • We shall adopt the following parameterization: 



TABLE I: Values of h(z) and e for which the studied metric 
is conformal to the indicated Bianchi types 



a 1 


h(6) = 1 - ez 2 (8) 


e 


Bianchi type 


dX 


sin 2 6 


1 


K-S 


dX 


1 





G 3 I 


dX 


sinh 2 9 


-1 


G3III 


dX + 21 cos 8dy 


sm 2 9 


1 


G3IX 


dX + 2lzdy 


1 





G3II 


dX + 21 cosh 6dy 


sinh 2 9 


-1 


G3VIII 



P12 = iH, P34 = D, 



(5) 



The studied line element is of the form 



where D is the electric displacement and H is the mag- 
netic field, the invariants acquire the form: 



f>ab p 
1 1 a., 



D 2 -H 2 



-, Q 



r>ab p 



iHD, (6) 



In this parameterization the structural function turns 
out to be 



dt 2 dz 2 
(j)(t) 2 ds 2 = — — + s(t)(dx + 2lzdy) + r—- + h(z)dy 2 , 
s(t) h(z) 

(10) 

where x and y are ignorable coordinates, h = 1 — ez 2 , with 
e and I being constants; s — s(t) and <j) — 4>(t) are the 
metric functions that are determined from the Einstein- 
Born-Infeld coupled equations. Chosing the null tetrad 
in the form: 



H 



^b i + b 2 (D 2 -H 2 )-H 2 D 2 
\J (b 2 + D 2 )(b 2 — H 2 ). 



(7) 



The parameterization JSJ has the advantage that there 
are only two nonvanishing components of the energy mo- 
mentum tensor: 



P 



:S4 



2b 2 



2b 2 



b 2 + D 2 



y/(b 2 + D 2 ){b 2 - H 2 ) 

b 2 - H 2 



sj(b 2 + D 2 ){b 2 - H 2 ) 
where we have defined: 



b 2 -H 2 
b 2 + D 2 



1] = 26 2 (e— - 1), 
1] = 2b 2 {e v - 1),(8) 

(9) 



The condition that e v be real imposes the restriction 
that H < b, i. e. the fields do not reach the maximum 
allowed field b. 



III. SPATIALLY HOMOGENEOUS 
SPACETIMES WITH NON-LINEAR 
ELECTROMAGNETIC FIELD 

Since the energy momentum tensor has two nonvanish- 
ing eigenvalues, Eqs. JHJ, may be coupled to spacetimes 
with two preferred directions, i. e. spacetimes of type D 
in the classification by Petrov. We have found solutions 
to the coupled Einstein-Born-Infeld equations for spa- 
tially homogeneous and anisotropic spaces, among them 
there are several cases of Bianchi metrics. 



i 



V2</> 
dt 



—= + iVhdy 



^2$ 
1 



—= + \fs{dx + 2lzdy) 



-—p + \/s(dx + 2lzdy) 
\/s 



The metric function s(t) turns out to be 



s(t) = + 2b l 



a -en 



dt 



dt 



(11) 



(12) 



where dot simbolizes derivative with respect to time i, 
C\ is an integration constant related to the vacuum case, 
b is the BI parameter and e is the constant curvature 
parameter. Moreover, <f>(t) must satisfy the differential 
equation <j) + l 2 <\> = 0. 

The system of EBI equations admits two possible so- 
lutions for the electromagnetic function v(t), Eq. 



2"W = y/l±<p*(t). 



(13) 



Both solutions present regions in which s(t) > pass- 
ing then to the opposite sign, s(t) < 0. The regions 
with s(t) > admit a cosmological interpretation as spa- 
tially homogeneous spaces. For s(t) < the signature of 
(fT7Tf> changes corresponding to stationary spaces (t di- 
rection changes from timelike to spacelike). We shall 
study the case e v = \J\ — 4 . The condition e u to be 
real imposes the restriction that (f> A < 1 and it determines 



4 



the range of the coordinate t. The electromagnetic fields 
are homogeneous and constant on spacelike hypersurfaces 
{t =const}: 



D(t) = (f> 2 (t) cos [21 



dt 



VM 1 

B(t) = -cp 2 {t)^[2l f-^=], (W) 
or in terms of the electromagnetic tensor: 

4?{t)F zy =B{t) = e v H{t), 



4> 2 {t)F tx =E(t) 



'D(jt), 



(15) 



The results in subsections III.B and III.C were derived 
in but are included here for completeness. 

The line element (11 (It can be written in terms of a 
synchronous time coordinate if we transform jkr 
and with a 1 = dx + 2lzdy: 



= dT 2 



dS 2 



-dT 2 + s(a 1 ) 2 + 



dz 2 



hdy 2 , 



(16) 



Written in the form (|16|) the Bianchi types can be iden- 
tified according to Table [Q(cf. Eqs. (13.1)-(13.2) in 19]). 
In case Bianchi G3IX the metric ()16[) is a particular case 
of Taub's universe (cf. Eqs. (6.29) or (8.5) in HJ). 



A. Bianchi cases included 

The spacetime ijTUjl possesses four Killing vectors. The 
four generator group G4 has invariant subgroups G3 that 
can be spatial rotations or translations (S3 or T3). In 
particular, we shall pay attention to the hypersurface- 
homogeneous spacetimes that can be identified as Bianchi 
spaces I, II, III, VIII and IX and one of them as 
Kantowski-Sachs. 

The Killing vectors X a are classified according to the 
values of the curvature parameter, e that can be 1, 0, — 1: 

i) e = 0, h = 1 



X\ — d y , X2 — d x , 



ii)e 



X 3 = -d z /2l + yd x , 

X A = zd y - yd z + l(y 2 - z 2 )d x 

-1, h = 1 + z 2 , 



(17) 



iii) e = 1, h = 1 — z 2 , 

X\ — d y , X2 = d x , 
X 3 — \fh cos yd z — SVCL ]l 



{2lyd x - zdy 



X 4 =Vhsinyd z + ^-(2lyd x - zdy) (19) 

Special cases of Bianchi metrics with a particular 
topology can be obtained as follows: 
i)e = 

In the metric Hl()[) with e = 0we obtain h(z) — 1 and 
the line element is 



dt 2 

cj>(t) 2 ds 2 = — + s(t)(dx + 2lzdyf + dz 2 + dy 2 , (20) 

s[t) 

where the ranges of all the coordinates are from —00 
to 00, this is a non-compact spacetime with topology 
R 4 that corresponds to Bianchi II. If I = 0, the metric 
becomes diagonal and can be identified as a Bianchi I, 



dt 2 

${t) 2 ds 2 = — ^ + s{t)dx 2 + dz 2 + dy 2 . (21) 
s(t) 

ii)e = -1 

In the metric (|10|l we transform z — > sinh# obtaining 
h{9) = cosh 2 6 and the line element becomes 



cj){t) 2 ds 2 



dt 2 
W) 



s(t)(dx + 2/sinh 9dy) 



dO 2 + cosh 2 9dy 2 



(22) 



where the ranges of all the coordinates are from — 00 to 
00, this is a non-compact spacetime with topology i? 4 
that corresponds to Bianchi VIII. If I = in the metric 
(|10fl and we transform z — > cosh6*, h(9) = sinh 2 9, then 
the line element is 



(j>{tfds 2 



d^_ 

s(t) 



+ s{t)dx 2 + d9 2 + sinh 2 9dy 2 , (23) 



that can be identified as a Bianchi III space. 
iii)e = 1 

In the metric IjlUII we transform z — > sa\9 and we ob- 
tain h(9) = cos 2 9 and the line element as 



X x 



d y , X2 — d x , 



X 3 = y/h cosh yd z - ^^^-(2lyd x + zdy), 

V 4 =Vhsivhyd z - ^^-{2lyd x + zd y ), (18) 



(j>(t) 2 ds 2 



dt 2 



+s(t)(dx + 2lsin9dy) 2 +d9 2 +cos 2 9dy 2 , 

where the ranges of the coordinates are — 00 < t < 00, 
< 9 < tt; if < y < 2tt, < x < 4?r it is a com- 
pact spacetime with topology Rx S 3 that corresponds to 
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Bianchi IX. For I = in metric (|10fl and changing coordi- 
nates to z — * cos 8 and with e = 1 we obtain h(9) = sin 2 9 
and the line element transforms into 



(j>{t) 2 ds 2 = + s(t)dx 2 + dO 2 + sin 2 9dy 2 , (25) 

s(t) 

it is a spatially-homogeneous solution which admits 
no simply-transitive G3 that can be identified as a 
Kantowski-Sachs space. 

B. Energy conditions 

According to Raychaudhuri, the equation that goberns 
the volume expansion of geodetic congruences, 9, is 



^ = -R ab V a V» + 2u? -2a 2 - U - + V a a (26) 
d\ 3 

where V a is a tangent vector to the geodesies, R a b is the 
Ricci tensor, w is the vorticity, a is the shear and V a is 
the acceleration of geodesies (or flow lines in a fluid); we 
shall consider that u = 0. 

From Eq. (12(it it can be seen that the expansion 9 
of a timelike geodesic congruence with zero vorticity will 
monotonically decrease along a geodesic if, for any time- 
like vector V a , R a bV a V b > 0. This is the strong energy 
condition (SEC), as stated in it says: "A spacetime 
satisfies the strong energy condition if R ab V a V b > 0, 
for all causal vectors V a ". The term involving the Ricci 
tensor, R ab , in Eq. lj2*6"|) induces contraction of the 
geodesic lines, indicating that the focusing of neighbour- 
ing geodesies requires SEC. 

For our Bianchi spaces the two nonvanishing compo- 
nents of the traceless Ricci tensor, in terms of the matter 
content, are 



Eqs. HI-©- DEC means that to any observer the local 
energy density appears nonnegative and the local energy 
flow vector is non-spacelike. In the BI case DEC is a re- 
quirement of a physically reasonable nonlinear structural 
function TL and its fulfilment is guaranteed by conditions 
(iii) and (iv) as stated in Sec. II. A. 

C. Scalar Polynomial Invariants 

For these solutions the only nonvanishing Weyl scalar 
is given by 




The scalar polynomial invariants constructed from \&2 
are given by / = S^/ 2 , and J = m such a man- 

ner that if ^2 is unbounded then both I and J must be 
unbounded and a curvature singularity occurs. The reg- 
ularity and boundedncss of ^2 depends on the proper for 
s\t),s{t) and s(t). In case of a suspected singular point, 
we must check if "J 2 is regular when evaluated at the sus- 
picious point. If ^2 remains bounded then it proceeds to 
investigate if the geodesies are complete. 

One could think that regularity of the curvature invari- 
ants might be helpful to avoid singularities but there are 
known examples of spacetimes with regular invariants, 
such as Taub-NUT, that enclose geodesies that are not 
complete in their affine parameterization due to the phe- 
nomenon called imprisoned incompleteness , in which a 
non-spacelike curve as it follows to its future enters and 
remains within a compact set. In spite that Hawking 
himself asserted that "the kind of behavior which occurs 
in Taub-NUT space cannot happen if there is some matter 
present", when studying a spacetime with regular invari- 
ants we must show that in fact geodesies are complete if 
one wants to assert that the space is singularity-free. 



Ri2 = 26 2 (1 - e"), i? 3 4 = 26 2 (l-e- y ). (27) 

If we consider the congruence V a = 4>^/sdt, or in the 
null tetrad formalism, V 3 — V± — A=, V 1 — V 2 = 0, for 

the BI field SEC amounts to R ab V a V b = i? 34 = 26 2 (1 - 
er v ) < 0, that clearly violates SEC. 

The violation of SEC is not exceptional, we remind 
that minimally coupled scalar field violates SEC and in- 
deed curvature-coupled scalar field theory also violates 
SEC. This point has been discussed recently in [21| and 
[22|. The aspect of interest for us is that the violation 
of this condition could imply the absence of an initial 
singularity since the fulfiling of SEC is related with the 
convergence of neighbouring geodesies. 

On the other side, the dominant energy condition 
(DEC) that states "for every timelike vector, V a , 
T a bV a V b > and T a bV a is a nonspacelike vector" 0, 
is indeed fulfilled by the BI energy-momentum tensor, 



D. Linear limit 

The linear (Maxwell) limit can be obtained from the 
solution Eq. (|12H taking the limit b — * 00, obtaining 

•v=*'^ +c2 Jw< + 'Jw ] - (29) 

where C 2 = H 2 + D 2 is related to the energy density of 
the field; in this case the electromagnetic field is given by 

H = Csin (2lat), D = Ccos (2lat), (30) 

where a is a constant. In Fig. ^ are displayed both met- 
ric functions, the linear (Maxwell) and Born-Infeld, the 
values of the constants correspond to a case of Bianchi 
G3 IX. In general the BI field smooths oscillations of 



6 



/ 

/ 
/ 

// 
if 

/ 
1 
\ 

// 


\ 

V ^ — 


' / 
/ 

V 





-1 1.5 3.5 



FIG. 1: The metric functions s(t) when the BI field is present 
(continuous curve) and in the linear limit (Maxwell field) 
(dashed curve) are shown. In general, the presence of BI 
field smooths oscillations of the metric function. The plot 
corresponds to Bianchi G3 IX, the values of the constants are 
I = 1 A = 0.8, C\ = 0.5, b = 0.85, e = 1 and the value of the 
constant linear field is C = 0.85. 

the metric function s(t). It is surprising that with lin- 
ear electromagnetic field as matter the spacetimes are 
singularity-free. We shall address geodesic completeness 
for this case in Sec. V. 



IV. ANALYSIS OF SINGULARITY 
STRUCTURE 

Hawking and Ellis Q discussed Bianchi I spaces and 
gave a theorem asserting that singularities will occur in 
all non-empty spatially homogeneous models in which 
the timelike convergence condition (R a f,u a u b > 0) is sat- 
isfied. However, as we showed in section III.B, in the 
studied case the timelike convergence condition is not 
fulfilled. Therefore the possibility of the absence of sin- 
gularity persists. 

For our Bianchi spaces with nonlinear electrodynamics 
there are two branches characterized by the value of the 
constant I that appears in the non-diagonal term of the 
metric, (dx + 2lzdy). In the case I = the conformal 
function <j)(t) fulfills the equation = 0, then <f> = At + B 



where A and B are constants. In the case I ^ the 
equation for <fi(t) is <f> + l 2 (f> = 0, with solutions <fi — 
A cos (It + B), A and B being constants. The analytical 
expression of the metric function s(t) depends on the 
conformal function <f>(t) and its behavior is different in 
each case. 

We will consider a spacetime as singularity-free when 
it is geodesically causal complete (g-completeness), that 
means intuitively that an observer in free falling does not 
leave the spacetime in a finite proper time; equivalently, 
that every geodesic can be extended to arbitrary values 
of its affine parameter. The demonstration of geodesic 
completeness is not simple in most cases Neverthe- 
less, the analysis of the geodesic extension can be simpli- 
fied when there exist constants of motion associated to 
Killing vectors. 

The method that we use to analize the completeness is 
to obtain the first order equations for the derivatives of 
the coordinates with respect to an affine parameter, x >Tl 
t,r, U,ti z .t, using the constants of motion. Then showing 
that these first derivatives are bounded, this implies that 
the corresponding geodesic curves x(t), t(r), y(r), z(t), 
are complete |23| . 

We shall analyze geodesic completeness of the space- 
times in the following subsections, first for the class of 
non-diagonal metrics, I ^ that includes Bianchi G3II, 
G3VIII and G3IX and then in the next subsection for 
the diagonal metrics Kantwoski-Sachs and Bianchi spaces 
G 3 I and G3III. 



A. Bianchi G 3 II, G3VIII and G 3 IX 

In the case that I 7^ the equation for <fi is <fi + l 2 <fi = 0, 
with solution tj> = Acos(lt + B), A,B being constants; 
the sector (t, x) of the metric is non diagonal and the 
cases included are Bianchi G3VIII and G3IX for e = 
— 1, 1, respectively and G3II when e = 0. 

The metric function s(t), Eq. I|12(l can be integrated in 
exact form; with <fr = A cos t (taking I = 1 and B = 0) we 
obtain a cumbersome expression in terms of sint, cost 
and elliptic functions, 



s(t) — — A^Ci cos 3 1 shit + e(cos t — cos 2 1 sin 2 1) 
2b^ 



, „ {cos 4 1 cos 2 t sin 2 t sin 2 t (2 cos 2t + cos it) J {2 - A 2 - A 2 cos 2 t)(2 + A 2 + A 2 cos 2 t) + 

A 2 3 3 6 



4- A 4 



8,2 4.2 (2 + A 2 + A 2 cos 2 1) 2 . o . , — r ^ 



F[7, /?]]}, 



. / v^sini \ . 2A 2 . , (I- A 2 )V2 + A 2 + A 2 cos2t\ 

a = arcsm — , a = — , 7 = arcsmh . (31 

\\J2 - A 2 - A 2 cos 2t J 1 + A 2 ' 1 \ 2sint )' V ' 



E and F are the elliptic integrals of second and first kind, respectively. From the expression l|31() . it is clear 
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x iT , since there will be divergences if s/<j) 2 is null, 
shall analyze this term with detail. 
Taking the limit t — > ir/2, 



lim 



We 



(34) 
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therefore, only in the case that e = 0, the geodesies 
x(t) can be incomplete. This case corresponds to spaces 
Bianchi II. Even in this case we can have complete 
geodesies when the constant K\ — 0, that means, 



FIG. 2: This is the characteristic plot for the metric func- 
tion s(t) and its derivatives, s(t),s(t) in the case I ^ 0. 
The three functions are periodic and bounded. The corre- 
sponding Weyl scalar, $2 does not diverge in the full range 
—00 < t < 00. In this graphic the values of the constants are 
A = 0.8, Ci = 0.5,6 = 0.85 and e = 1. For e = 0, -1 the 
behavior is qualitatively the same. 



that s(t) is a periodic and bounded function. There are 
values of t for which the metric function s{i) is null and 
there the metric (I1U|I is singular. For <f>(t) = ^4cos(i), 
these values are t = ±(2n + l)ir/2,n = 0,1,2,... and 
there we might have a singularity. There are additional 
zeroes when negative terms cancel the positive ones, how- 
ever, the position of these zeroes depends on the balance 
between the distinct parameters, Ci,b 2 /A 2 , A 2 and e. 

The characteristic behavior of the functions 
s(i), s(t), s(t) is illustrated in Fig|3 The three functions 
s(t), s(t), s(t) are periodic and bounded, that do not di- 
verge over all the range —00 < t < 00. As a consequence, 
the Weyl scalar ^2, Eq. (|28|l . is also finite over all the 
range — 00 < t < 00 and so are its derived invariants. To 
dilucidate the behavor of the spacetime at t = tt/2 we 
shall invetigate the completeness of geodesies. 

For each Killing vector Xj there exists a conserved 



Ki = -^(x,r + My, T ) =0, 



(35) 



quantity x^X^ 



v^.g^X" =const, where x m (t) is a 



geodesic parametrized by r. Since d x and d y are Killing 
vectors, it implies that there are constants K\^K% given 
by 



K\ = g xx x tT + g xy y,T, 
K 2 = g yy y,T + g yx x,T, 

From Eqs. (|32|) we determine x tT and yy. 



-{K 2 -MK{), 



(32) 



My-, 



(33) 



where M — 21 z. From Eqs. (|33|) we can analyze com- 
pleteness for x(t) and y(r). All the coefficients in Eqs. 
(|33|l are finite (h = 1 — ez) except for the first term in 



in that case the geodesies x. T = —My tT) but the equation 
for y >T is given in terms of bounded coefficients, also M = 
21 z is bounded, then in this particular case of e = 0, 
the geodesies x(t) are complete. However the generic 
behavior when e = is geodesically uncomplete. 

To obtain the equations for t_ T and z iT we shall use the 
rest of the Killing vectors, quoted in Sec. III. A. These 
Killing vectors are different depending on the value of 
the curvature parameter e. We separate in three cases: 

6=1,-1,0. 

For e = 1 the Kiling vectors Eqs. (|19fl imply the exis- 
tence of two conserved quantities, K3 and K4 that permit 
determine z T as 



z 2 T + (f) 4 (zK 2 - 21K 1 ) 2 - <f) 4 h{K 2 + K 2 ) = 0, (36) 

On the other hand, we can also use the line element, 
to obtain i jT , 



ds 2 



r^-:£ + -A^ + A Ar) 2 , (37) 



where S is 1 for spacelike geodesies, for null and -1 if 
the geodesies are timelike. 

Substituting the constants of motion we have, 

t 2 T = s^{(K 2 + Kl)-(zK 2 -2lK 1 ) 2 /h + 

(K 2 ~ KxMf/h} + <j) 4 Kf - <f> 2 sS. (38) 

In the equations for t 2 T and z 2 T all the coefficients are 
bounded, then the geodesic curves t(r) and z(t) are com- 
plete. 

In an analogous way can be determined t }T and z T for 
the case e = — 1. Using the Killing vectors X 3 and X4, 
Eqs. fl%|l we obtain: 



(zK 2 + 21K X ) 2 + cj> A h{Kl - K : 



t 2 T = sct) i {(Kl-Kl) + (zK 2 + 2lK 1 ) 2 /h + 
{K 2 - K^lf/h} + <fK\ - <j) 2 s8. 



(39) 
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FIG. 3: It is shown the plot for the metric function s(t) 
and its derivatives, s(t),s(t) for the case I = 0. The three 
functions are bounded at t = 0, but s(f)diverges at at — > 1. 
In this graphic the values of the constants are A = 0.8, C\ = 
0.5,6 = 0.8 and e = 1. For e = 0, —1 the behavior changes 
quantitatively. 



For e = we have from the motion constants that 
imply the Killing vectors in Eqs. (|17fl . 



Al(j) 4 h 2 {lK 2 - K4K1 - KxilKxz 2 - K 2 z)}, 
S( f> 4 {Al 2 h{K 3 - K lV ) 2 + [K 2 - KtMf/h} 



l K 2 - tfsS. 



(40) 



In the three cases, e = 1, — 1, the geodesies t(r) and 
z(t) are complete. 

If one considers the topology RXS 3 for Bianchi G3IX 
(case e = 1) then < x < Air and x reaches a finite 
value as r increases. This geodesic "wraps" around the 
x-direction an infinite number of times. It corresponds to 
the "phenomena" known as imprisoned incompleteness, 
described by Hawking (see paragraph 8.5 in P) for Taub- 
NUT spaces. This pathological behavior arises essentially 
because x is identified at and 4-7T, a consequence of the 
compactness of Bianchi type IX. 

However both types IX and VIII can be considered 
non-compact and we can take the range of x being 
(—00, 00). Being x no longer periodic the geodesies x(t) 
can be extended indefinitely and the spaces are complete. 
The space VIII resembles the ones studied by Siklos [24| 
as extensions of Taub-NUT metrics. 



s(t) 



Ab 2 t 3 {K[ 



2V2 
Oh 2 , 



— arcsin At, -]} 
dA^-et 2 , 



(41) 



where F is the elliptic integral of the first kind and K is 
the complete elliptic integral of the first kind. We observe 
from Eq. 141|) that s(t) take real values in the interval 
< (At) 4 < 1 and that s(t = 0) = 0. 

The metric function s(t), and its derivatives, s(t) and 
s(t) are finite at t — 0, then the Weyl scalar Eq. lt2S|l . 
and the invariants derived from it are regular at t = 0. 
This is an indication that singularity at t = is not a 
physical one. However, s(t) diverges as At — * 1 and it 
implies that the polynomial invariants constructed from 
^2 diverge as well at At = 1, having there the possibility 
of a singularity. We shall analyze the geodesies to figure 
out the singularity structure if it exists. Fig. displays 
the metric function s(t), and its derivatives, s(t) and s(t). 

The equations for the first derivatives of the coordi- 
nates with respect to the affine parameter r are equally 
valid in this case, just putting I = 0, M — and the met- 
ric function s(t) from Eq. (|41|) . If the first derivatives 
are bounded it implies that geodesies are complete and 
the field is nonsingular [2]|; the first equations being, 



y.r 



h : 

K^ 2 



(42) 

The term 2_ that might diverge, is bounded in all the 
range of t and diverges at t = only when e = 0, because 



lim — — = — 



c 

i 3 ' 



(43) 



in such a manner that geodesies x(t) are unbounded for 
Bianchi G3I (e = 0), since x^ T diverges for a finite value 
of the affine parameter. Even for the case e = if the 
constant K\ = then the geodesies x(t) are complete 
since in this case x(t) = At + B, A, B constants, and it 
is a complete curve. 

The coefficients in equations for z r and t tT are 
bounded as can be seen from expressions H36JH4()|I 
putting / = Q, M = . On the other side, in spite 
that 



B. Kantowski- Sachs and Bianchi spaces G3J and 
G3III 

Let us address the case 1 = 0. The included spaces are 
Kantowski-Sachs and Bianchi spaces G3I and G3III for 
e = 1,0,-1, respectively. 

In this case (f> = At is only a scaling of time coordinate; 
s(t), Ea. l|12|l . can be integrated in exact form as 



lim (si 

ai->l v 



(44) 



all the coefficients in the equations for x a r , x a — (x, y, t, z) 

,2 

are bounded, in particular the term — is finite there, 



lim — = C X A 



2V2b 2 1 2b 2 e 
^ K[ 2 ] ~3A^-A^- (45) 
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Then at (at — > 1) the geodesies are complete, showing 
that there is no singularity at that time. 

In relation to the electromagnetic fields, it can be seen 
from Eqs. i|14fl that when I = only electric or only 
magnetic field persists (remind the electromagnetic du- 
ality, which is preserved in the BI theory), for instance, 
D(t) = (f>{t) 2 , B = 0. But note that the fields can 
not trespass the threshold of the maximum field b, then 
D(t) — A 2 t 2 < b imposes an aditional restriction for the 
time coordinate, it must be that At <Inf {Vb, 1}. 



V. SINGULARITY-FREE SOLUTIONS WITH 
LINEAR ELECTROMAGNETIC FIELD 

The conserved quantities derived from the existence of 
the four Killing vectors are independent of the linear or 
nonlinear electrodynamics. Then the equations for the 
derivatives of the coordinates with respect to the affinc 
parameter can be analyzed to test completeness if instead 
of having the BI field, we take the linear limit. 

Expressions for the derivatives of the coordinates with 
respect to the affine parameter, Eqs. (|23J)-(@UJ) are valid 
now with s(t) given by its linear limit, Eq. Ij29|l . The 
possibility of divergence appears again in the term s/(j> 2 
in the equation for x tT . We analyze the behavior of the 
term s/<j> 2 . Integrating the expression for s(t) Eq. i|29|) 
with d> = A cos t we obtain: 
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— [ 

A 4 



s (t) = C 1 ^+^ 

(46) 

where C 2 = H 2 +D 2 is proportional to the energy density 
of the field and the limits at t — > and t — > ir/2 are 



lims(t) = C l IA 2 + e, 



lim s(t) 



-C 2 /3A 2 



(47) 



From the previous expressions it is clear that s(t) is 
not zero neither at t = nor t = 7r/2, being the pres- 
ence of the electromagnetic field the term that allows the 
completeness of the geodesies; therefore the coefficients 
in the equation for bounded and consequently the 

geodesic x(t) is complete. This proves that for the solu- 
tion of the Einstein-Maxwell equations given by the line 
element (|1U|I with s(t) in Eq. I|4()|) there are singularity- 
free spacetimes. 

Moreover, for the case I = we also have singularity- 
free solutions since the term s/4> 2 is not zero neither at 
t = not at At = 1; in this case <j) = At and 



s(t) = dAH 3 -C 2 -et 2 , 
and the limits 



(48) 



lims(t) = -C 2 , 

lim s(t) = CiA-C 2 - ej A 2 , 

at—>-l 



(49) 



Showing that the geodesies are complete, since the co- 
efficients in the equation for X t 3X6 bounded. Summa- 
rizing, there are singularity-free solutions with constant 
electromagnetic field for spacetimes given by the metric 
(|10fl and the metric function s(t) given by (|46[) if I ^ 
or by H48|) if / = 0. The electromagnetic field is given by 
Eqs. i|3U|) . We also note that for these singularity- free so- 
lutions the electromagnetic field tensor does not violate 
SEC. 



VI. CONCLUSIONS 

We investigated spacetimes with a four-dimensional 
group of isometries with a three-dimensional subgroup 
acting transitively on spacelike hypersurfaces, i. e. spa- 
tially homogeneous spacetimes anisotropic in one spatial 
direction. We address the study of geodesic completeness 
for solutions to nonlinear electromagnetic field coupled 
with Einstein equations for cases of Bianchi G 3 I, G3II, 
G 3 III G 3 VIII, G 3 IX and Kantowski-Sachs. 

Two families are distinguished: Bianchi G 3 I, G 3 III 
and Kantowski-Sachs on one side and Bianchi G 3 II, 
G 3 VIII and G 3 IX on the other. 

Bianchi G 3 I, G 3 III and Kantowski-Sachs only admit 
one component of the electromagnetic field, electric or 
magnetic (not both, interchanging roles by a duality ro- 
tation). These spaces show divergence in their invariants 
at a finite time, however, geodesies are complete. At 
t = it is shown that geodesic uncompleteness appears 
for the case e = 0. 

Types VIII and IX are geodesically complete and as 
such they are singularity-free. The space VIII resembles 
the ones studied by Siklos [24( as extensions of Taub-NUT 
metrics. When Bianchi G3IX is considered as having the 
topology RXS 3 it presents imprisoned incompleteness, 
described by Hawking (see paragraph 8.5 in 0) for Taub- 
NUT spaces. This pathological behavior arises essentially 
as a consequence of assuming compactness for Bianchi 
type IX. 

Our results are not in contradiction with the assertion 
by Hawking about the necessity of singularity in non- 
empty spatially homogeneous models, since the timelike 
convergence condition fails to hold in nonlinear electro- 
dynamics. The incompleteness due to the presence of 
imprisoned geodesies in Bianchi G3IX seems to be incon- 
sistent with the affirmation that: in the Taub-NUT case 
the incompleteness of geodesies leads to a real singular- 
ity once any amount of matter is added (see for instance 
[23). In this paper we have presented a counterexample 
for a space of Taub-NUT type with imprisoned geodesies 
but there is no contradiction since we have considered as 
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matter a nonlinear electromagnetic field and the original 
assertion was made thinking of matter as a perfect fluid. 

Finally we remark that in the linear limit (Einstein- 
Maxwell) the solutions are singularity-free for all the 
Bianchi cases included, being clearly in this case the pres- 
ence of the field that makes the geodesies complete. In 
the nonlinear case the determinant factor to avoid the 
singularity is the curvature parameter, existing singular 



[1] Hawking S and Ellis G 1973 The large scale structure of 

space-time, (Cambridge University Press, Cambridge) 
[2] Geroch R 1968 Ann. Phys. 48, 526. 

[3] Mars M and Senovilla J M M 1997 Class. Quantum Grav. 
14 205. 

[4] Berger B K 2000 Phys. Rev. D 61 023508. 

[5] Misner C W 1969 Phys. Rev. Lett. 22 1071. 

[6] Toporensky A V and Ustiansky V O Dynamics of 

Bianchi IX universe with massive scalar field: Preprint 

gr-qc/9907047 

[7] Garcfa-Salcedo R and Breton N 2000 Int. J. Mod. Phys. 
A 15 4341. 

[8] Garcfa-Salcedo R and Breton N 2003 Class. Quantum 
Grav. 20 5425. 

[9] Jackson J D 1975 Classical Electrodynamics Second Ed. 
(John Wiley and Sons). 
[10] Heisenberg W and Euler H 1936 Z. Phys. 98 714. 
[11] Drummond I T and Hathrell S J 1980 Phys. Rev. D 22 
343. 

[12] Schwinger J 1951 Phys. Rev. 82, 664. 

[13] A. A. Tseytlin A A 1996 Nucl. Phys. B 469 51; Fradkin 



solutions if the curvature parameter is zero (e = 0). 



Acknowledgments 

R. G-S. acknowledges partial support from project 
UAEM 1941/2004. 



E S and Tseytlin A A 1985 Phys. Lett. B 163 123. 

[14] Brill D R 1964 Phys. Rev. 133 845. 

[15] Hughston L and Jacobs K 1970 Astrophys. J. 160 147. 

[16] Hacyan S 1984 Astron. Astrophys. 140 363. 

[17] Plebanski J F 1970 Lectures on Nonlinear Electrodynam- 
ics (Copenhagen). 

[18] Plebanski J F 1973 Spinors, Tetrads and Forms (Ed. 
CINVESTAV-IPN, Mexico). 

[19] Stephani H, Kramer D, MacCallum M, Hoenselaers C, 
Herlt E, 2003 Exacts Solutions of Einstein's Field Equa- 
tions. (Cambridge University Press). 

[20] Ryan M P and Shepley L C 1975 Homogeneous Relativis- 
tic Cosmologies (Princeton Univ. Press, Princeton). 

[21] Senovilla J M M 1988 Gen. Rel. and Grav. 30 5701. 

[22] Visser M and Barcelo C 2000 Energy conditions and their 
cosmological implications Preprint gr-qc/0001099 

[23] Arnold V I 1990 Ordinary Differential Equations, (MIT 
Press) . 

[24] Siklos S T C 1976 Phys. Lett. A 59 173. 



